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The measurement of quarkonia production in relativistic heavy ion col¬ 
lisions provides a powerful tool for studying the properties of the hot and 
dense matter created in these collisions. To be really useful, however, such 
measurements must cover a wide range of quarkonia states and colliding 
species. The PHENIX experiment at RHIC has successfully measured 
J,Xc and Upsilon production in different colliding sysyems at vari¬ 
ous energies. In this talk I will present recent results from the PHENIX 
collaboration on charmonium production in d-|-Au, Cu-|-Au and U-l-U col¬ 
lisions at 200 GeV/c. 
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1 Introduction 


Dissociation of quarkonia by color screening in deconfined matter is predicted to be 
different for different states. Loosely bound states will melt first, and successive 
suppression of individual states can provide an effective thermometer of the QGP. 
However, there are many competing processes in nucleus-nucleus collisions: cold nu¬ 
clear matter effects, color screening, initial state effects, regeneration, feed-down, and 
so on. Thus, in order to have a clear picture of what happens during relatistic nucleus- 
nucleus collisions, we need measurements for different energies, colliding species, and 
quarkonium states. In this respect, asymmetric nucleus-nucleus collisions can be very 
useful in understanding the importance of different processes contribuing to quarkonia 
production. 

2 d+Au collisions 

The PHENIX experiment has measured J/ip production in p+p and d-l-Au collisions 
at 200GeV. at forward, backward, and central rapidities jl]. Fig. [I](top left) shows 
the J/?/^ invariant yields in p-|-p and d-l-Au collisions as a function of rapidity, in¬ 
tegrated over centrality (0%100%). The error bars (boxes) represent point-to-point 
uncorrelated (correlated) uncertainties. 

The cold nuclear matter effects are quantified by calculating the nuclear modifi¬ 
cation factor RdAu) which is defined as the ratio oi J/'ip yield in d-l-Au collisions to 
J/'i/i yield in p-|-p collisions, corrected for the number of binary collisions Ncou- Ncoii 
is derived using a Glauber calculation (see [8] for details). 

As expected, nuclear modification factor RdAu exibits rapidity asymmetry in d-l-Au 
collisions, as is shown in Fig.[I](bottom left). Forward (deuteron going) rapidity shows 
more suppression than central and backward (Au going) rapidity. 

This difference can have many possible explanations, including nuclear breakup 
and gluon shadowing. A model [6] which uses EPS09 nPDF and breakup cross-secton 
abr — ^mb shows reasonable agreement with the data. Predictions of this model are 
shown in Fig. [T] by red curves. A second class of calculations incorporates gluon 
saturation effects at small-x [7j, and is shown in Fig. [D by green lines. This model 
the data well at forward rapidity, but fails to reproduce them at backward rapidity. 

However, if one looks at centrality dependence of RdAu the picture becomes less 
unambiguous. Fig. [fright) shows rapidity dependence of RdAu for peripheral (top) 
and central (middle) d-l-Au collisions, as well as their ratio Rqp (bottom). For pe¬ 
ripheral collisions, the RdAu ratio shows a mild suppression, roughly independent of 
rapidity, within the systematic uncertainties of approximately 15%. For central col¬ 
lisions RdAu indicates a much larger suppression for J/tp at forward rapidity. As one 
can see, the model which uses EPS09 nPDF and breakup cross-secton abr = 4m6 [6] , 
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Figure 1: Rapidity depencence oi J/'ip yield in p+p and d+Au collisions (top left), 
and of RdAu (bottom left). Rapidity dependence of RdAu for different centralities 
(right). 


fails to describe the Rcp measurement at forward rapidity. No parameter choice of 
the EPS09 nPDF set and of abr is able to describe the rapidity and centrality depen¬ 
dence of the data. Gluon saturation model [7] describes centrality dependence well 
at forward rapidity, but fails at other rapidities. 

Fig. [2] shows transverse momentum dependence of RdAu for different rapidities. 
As one can see, at all rapidities, RdAu rises up to 5GeV/c. A model which includes 
shadowing -|- abr, shown by red dashed curves, does not match the trend. The model 
by Kopeliovoich et al. [S], which includes Cronin effect and abr qualitatively matches 
the shape of observed dependence. Largest disagreement with theories is observed at 
backward rapidity. 
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Figure 2: Transverse momentum dependence of RdAu for three different rapidity 
ranges. 
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Fig. El shows RdAu for [3]. Unexpectedly, is 3 times more suppressed in 
most central collisions than J/ip, and has very different trend with N^oii compared to 
J Reference [9] presents a model that explains the lower energy E866/NuSea and 
NA50 relative suppression results using an expanding color neutral cc pair. 

As the cc expands, it has an increased nuclear absorption owing to its larger physical 
size. Once the time spent by the cc pair traversing the nucleus (r) becomes larger 
than the Jxp formation time, the ip' will see a larger nuclear absorption owing to its 
larger size. This idea is tested at RHIC energies by calculating the average proper 
time spent in the nucleus by the quarkonia (or cc precursor). 

Fig. m shows nuclear crossing time in d+Au for different collision energies. Uni¬ 
versal trend with dNch/drj for several systems, up to 200 GeV is observed. The solid 
curve in Fig. E] is the calculation by Arleo et ah jH], which is consistent with the trends 
observed by E866/NuSea and NA50. However, the PHENIX data show very different 
T dependence. 

The values of r for the PHENIX data are similar to the cc formation and color 
neutralization time of ~ 0.05/m/c, and well below the Jip formation time of ~ 
0.15/m/c [H]. Therefore the model cannot explain the strong differential suppression 
of the tp' in the PHENIX data. 
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Figure 3: RdAu for ip' and J/ip as a function of centrality. Note that the J/ip RdAu 
plotted here is not corrected for Xc and ip' feed-down, and the Ncoll values are shifted 
slightly to aid in clarity. 


3 Cu+Au collisions 

Fig. E] shows J/ip Raa vs. centrality for Cu-|-Au collisions [TU] as open black cir¬ 
cles (backward rapidity) and solid black circles (forward rapidity). Raa for Au-|-Au 
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Figure 4: Relative modification ^’/ J/^ as a function of quarkonia time in nucleus. 
The curve is a calculation by Arleo et al. [H] discussed in the text. 

collisions is shown on the same plot for comparison as orange points. Observed sup¬ 
pression is somewhat smaller in peripheral Cu-|-Au collisions, compared to Au-|-Au 
collisiosn, but becomes the same with increasing centrality. Higher suppression is 
observed in the region of lower particle density (forward rapidity), similar to d-l-Au 
collisions. Debye screening would go in the other direction. 

Fig. El shows a ratio of Raa for Cu-going and Au-going directions. This ratio has 
the advantage of reduced systematic uncertainties. Observed ratio decreases with 
centrality. The 20% 30% difference in suppression between forward and backward 
rapidity Raa evident in Fig. |6] could be due to hot matter effects, CNM effects, or 
a combination of both. Also shown in this figure as solid line is a model m which 
estimates the contribution from cold nuclear matter (shadowing). The grey band in 
this figure represents the extreme nPDF parameter sets for the model. 

As can be seen form Fig. [6] The difference between forward (Cn-going) and back¬ 
ward (An-going) J/^/j modification is found to be comparable in magnitude and of 
the same sign as the expected difference from shadowing effects. 

4 U-|-U collisions 

Collisions of deformed uranium nuclei produce a wide variation in energy density 
within the same colliding system. MC studies show [T] a possibility of selecting 
experimentally tip-tip collisions by selecting high multiplicity, but low flow events. 
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Figure 5: Raa vs. centrality for Cu+Au and Au+Au collisions. 



Figure 6: Raa ratio for Cu-going/An-going directions in Cn+Au collisions. 


In tip-tip collisions T/Tq could reach above 2 [2], at which tempertaure T(IS') could 
dissociate. 

The PHENIX collaboration recently measured J/ip production in U+U collisions 
at 200GeV. Fig. [3 shows preliminary results of J/ip Raa as a function of central¬ 
ity for U+U collisions (black circles) compared to Au+Au and Cu+Cu collisions. 
Qualitatively similar J/ip suppression is observed from Cu+Cu to U+U collisions. 
Somewhat weaker suppression in most central U+U collisions may indicale higher 
role of coalescence processes in nraninm-nraninm collisions. 
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Figure 7: Raa vs. centrality for U+U collisions (black circles) compared to Au+An 
and Cu+Cu collisions (open symbols). 

5 Concisions 

In d+An collisions J/%1) nnclear modification factors at forward rapidity as a fnnc- 
tion of centrality cannot be reconciled with a pictnre of cold nuclear matter effects 
(nPDFs and when an exponential or linear dependence on the nuclear thickness 
is employed. 

In Cu+Au collision, the Cu going side is more suppressed than Au going side, 
consistent with CNM effects (shadowing). 

The magnitude and trend of V'(2s) suppression in nnclear collisions is qnite dif¬ 
ferent from ithat of J/ip. Nuclear crossing time does not explain the data. 

J/ip Raa is qualitatively consistent between different colliding systems, from 
Cu-|-Cu to U-l-U. ~25% differences could be due to expected variations in the CNM 
and QGP effects. 
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